NASA-CR- 194648 


•// 


U 



/f/7i 

Final Report to the National Aeronautics and Space Administration 
for the ROSAT Grant entitled: 



THE FORMATION PROCESS OF THE HeT A10830 LINE IN 

COOL GIANT STARS 


NAME AND ADDRESS OF Iowa State University 

INSTITUTION: Ames, Iowa 50011 


NASA GRANT NUMBER: NAG 5-1655 

ISU ACCOUNT NUMBER: 422-21-04 


PRINCIPAL INVESTIGATOR: 


REPORT COVERING THE PERIOD: 


Dr. Donald G. Luttermoser 
Iowa State University 

7/01/91 - 6/30/93 


-4- 

o 

o 

r— i/i 

rt 

* o 

o c 

Z 3 


O 

a 

\ 

m 

O 


m 


o 

*-4 

o 


DISTRIBUTION: 


Three Copies to: 

Dr. Robert Petre — Code 666 
NASA ROSAT Technical Officer 
Lab. for High Energy Astrophysics 
NASA/ Goddard Space Flight Center 
Greenbelt, MD 20771 

Two Copies to: 

NASA Sci./Tech. Info. Facility 
Accessioning Department 
800 Elkridge Landing Road 
Linthicum Heights, MD 21090 

One Copy to: 

Mr. Richard Hasbrook 
Office of Contracts and Grants 
Iowa State University 
204 Beardshear 
Ames, Iowa 50011 


m 

o 

O Ch 

2 rH 

C? cc p— 

O rv • 

*-< c c 

< •• D 

X < LL “■> 

Q 

C a t/> o 

t-i- X XT 
< < 


m 


OJ 

X 


c 

q 

u 

c 

0) 

V 

CO 


CO 


I 


h— r — 4 O 

K O 

v # 

^ JL <T r-4 ^ 

t) WH ... 


a 


sr UJ (' * CT r—4 

OX *—3 

st I- J ^ 

O' > 1 / 

— • a o >■ 

• ..J o ro 

o- O 

0 co * ./> _ 

1 to C 

< 'Jj C nj c 

CO CO u O ’fcJ" 

< 3 Cl b U 

' CL X ' K* 


One Copy to: 

Ms. Gloria R. Blanchard, Grants Officer, 
NASA/ Goddard Space Flight Center 
Greenbelt, MD 20771 


REPORT SUBMITTED: 


October 29, 1993 
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The Formation Process of the He I A10830 Line in Cool Giants 
Donald G. Luttermoser, P.I. (Iowa State University) 

This NASA grant covers my research program, which involves observing a sample of cool giant 
stars with ROSAT. These stars were selected from the list of bright stars which display He I 
A10830 in absorption or emission and lie on the cool side of the coronal dividing line. With 
measured X-ray fluxes or upper limits measured by the PSPC, we investigated the role X-rays 
play in the formation of this important line using the non-LTE radiative transfer code PANDORA. 
Hydrodynamic calculations have been performed to investigate the contributions of acoustic wave 
heating in the formation of this line as well. 


1 Introduction 


Linsky & Haisch (1979) determined that the cool side of the H-R diagram is divided into 2 distinct 
regions: (1) F-M dwarf and late F through early K giants, which show evidence for chromospheric 
( T < 20, 000 K) and higher T (5 X 10 4 — 2 x 10 5 K) gas through emission lines seen in IUE spectra, 
and have been coined coronal stars ; and (2) Single red giants later than K2 III and supergiants 
later than G5 lb, which have generally weaker chromospheric emission and no evidence for higher 
T species — the so-called non-coronal stars. The “line” that separates these regions has been 
coined the Coronal Dividing Line (CDL). Ayres et al. (1981) confirmed the CDL with Einstein 
observations at X-ray wavelengths. Stars on the cool side of the CDL also show evidence for massive 
stellar winds (10 -8 — 1 O - 5 M 0 yr -1 , Stencel & Mullan 1980). These observations have lead to the 
suggestion that a “mechanical” energy input on the hot side of CDL heats gas to temperatures 
greater than 10 6 K, and that a mechanical energy flux on the cool side of CDL drives strong stellar 
winds. 

Coronal stars are thought to be solar-like with: ( 1 ) Chromospheric layers where 


Imsh. > o —>o 

dh ' dh ’ 


T ~ 10 4 K 


and 


h corona > h 


which give rise to emission features in singly ionized resonance lines (i.e., Ca II H & K, Mg II 
h & k, Fe II). The strength of ionized emission features have been noted to be correlated with 
5-field strength (see Schrijver 1991). This suggests that acoustic heating of sound waves from the 
convection zone of these stars is modified (somehow) by magnetic fields. ( 2 ) Transition Region 
(TR) layers where 

Fmech > 0 , ^ » 0 , T - 10 5 K 


dh 


and 


dh 


h corona ^ h > 


which give rise to emission features in doubly- and triply- ionized resonance lines (?.e., C IV A 1550). 
(3) Coronal layers where 

%^>0, f.o, T ~ 2 x 10 6 K 
dh dh 
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and represent the uppermost layer, which give rise to emission features of highly ionized species at 
X-ray wavelengths and show a definite correlation between X-ray flux and 5-field strength. 

Meanwhile, non-coronal stars have no evidence of magnetic field influence, but instead have ev- 
idence for: (1) Chromospheric layers which give rise to emission features in singly ionized resonance 
lines; and (2) a wind region. Much debate has ensued over the characteristics of the chromospheres 
in these late-type stars. Are they solar-like, a hydrodynamic chromosphere or “calorisphere” 
(Bowen 1988; Willson 1988; Dupree et al. 1990; Luttermoser & Bowen 1992), or possibly a post- 
shock region (Bowen 1988; Luttermoser & Bowen 1990)? There is also considerable discussion 
concerning the mechanism that drives the winds of these stars. Possibilities include radiation pres- 
sure on dust, molecular levitation (Elitzer et al. 1989), and shock propagation and merging (Bowen 
1988; Cuntz 1992). 

Stars on the “hot” side of the CDL display He I (2s 3 S - 2 p 3 P°) triplet lines near 10830 A that 
are broad, shallow, and constant, similar to j3 Gem (O’Brien & Lambert 1986). The formation 
of this chromospheric multiplet in these stars is controlled by backflowing coronal EUV and X- 
radiation. Lambert (1987) has shown that stars on the “cool” side of the CDL show variable 
He I A10830 features — sometimes going into emission! With this grant, we wished to answer the 
question: “How is this line formed in the non-coronal stars?” 

Einstein detected no X-rays from nearby K-giants Fx < 0.017^(0) (Ayres et al 1981). The 
ROSAT observations of this program and others have pushed this limit to Fx < O.OOlEv(O) for 
a Tau (Luttermoser 1994) and Fx < 0.00017^(0) for a Boo (Ayres et al 1991). Should this 
flux be too low to produce the observed He I A10830 feature via photoionizations, one would need 
T > 40, 000 K gas for collisional excitation to populate the 2s 3 S state as shown in Figure 1 (Cuntz 
& Luttermoser 1990). Stochastic shocks formed by the propagation of short period acoustic waves 
into the tenuous layers of the outer atmosphere are a likely scenario (Cuntz 1987). 

2 ROSAT Observations of a Tau & a Cas 

Two stars just on the “cool” side of the CDL (see Figure 2) were selected as a test case in under- 
standing He I A 10830 formation in these types of stars. The star a Cas is classified as K0 II— III 
(Le., a 4700 K star at a luminosity class between bright giant and giant). It has an apparent, 
integrated-bolometric flux hoi — 4.6 x 10“ 6 erg cm -2 s -1 and d — 63 parsecs. Meanwhile, a Tau 
is a K5 III (i.e., 4000 K giant) star with hoi = 3.6 x 10" 5 erg cm -2 s' 1 and d - 20 parsecs. These 
stars both show He I A10830 lines. This line in a Cas has W\ = 230 mA and is constant in strength 
similar to f3 Gem (Lambert 1987). In a Tau, it varies in strength from no line to W\ = 37 mA. 
It was not observed in emission however as is the case for some non-coronal K giants (O’Brien & 
Lambert 1986). 

IUE observations have suggested that both of these stars have chromospheres but no coronae 
(see Ayres et al 1981). For a Cas: fM g ri/hol = 8.2 X 10~ 6 with no apparent emission in C IV or 
He II (A1640). For a Tau: fM g ii/(boi — 4.8 x 10 -6 with no apparent emission in C IV or He II. For 
reference, the known coronal star (3 Get (K1 III) has f Mg iilhol = 2.6x10 ~ 5 ,fciv/hoi = 1.0 xlO" 7 , 
fHell/tbol = 1.1 x 10“ 7 . 

The Einstein observations by Ayres et al. (1981) did not detect any X-ray flux for a Cas or 
a Tau. Upper limit for a Cas were fx /hoi < 3.0 X 10 -8 and for a Tau were fx/hol < 4.0 X 10“ 9 . 
For reference, / 3 Get has fx/hol — 4.0 x 10 -6 . The K0 III star (3 Gem, for which the He I A10830 
line does not vary, has fx /hoi = 4.0 X 10 -8 , and represents the weakest coronal X-ray emitting 
giant star seen by Einstein . 

In order to understand the mechanism of the formation of the He I A10830 feature for stars on 
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Figure 1: Electronic structure of He I. Note that the lower state of the A10830 multiplet is a 
metastable state. 
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Figure 2: vs - V — R diagram for a sample of stars near the CDL. Our target stars a Cas and 

a Tau are indicated. Note that a Cas is on the “non-coronal" side of the CDL yet has a He I 
A10830 signature similar to a “coronal” star. 




the cool side of the CDL, we had to first push down the Einstein upper limits of the X-ray flux 
for these stars. During ROSAT’s 1st and 2nd years, the PSPC (Position Sensitive Proportional 
Counter) obtained data over 3 separate pointings for a Tau on 24 February 1991 (Figure 3) over a 
total of 6833 seconds, and a Cas on 7 January 1992 (Figure 4) over a total of 5897 seconds. The 
PSPC has an energy range of ^0.1-2. 4 keV (124-5 A) and a modest degree of energy resolution 
(AE/E ^ 0.5) (see Triimper 1983). The PSPC did not detect a Tau in 7 ksec (note that Tom 
Ayres did not detect a Boo (K1 III) in 18.6 ksec). 

Using IDL reduction software developed by Dr. Fred Walter, the new X-ray flux uj^per limit, 
for a Tau is 2.08 X 10" 15 erg cm -2 s" 1 , which is based on the 3 <r upper limit of 4.8 counts. 
This corresponds to fx/hol < 5.78 X 10~ u (a factor of 70 improvement to the Einstein upper 
limit). Previously undetected in X-rays, this program discovered a Cas to be a weak X-ray 
source! The net extracted source count is 56.6 ± 8.7 during this observation. This corresponds to 
fx — 2.46 X 10" 14 erg cm" 2 s" 1 (note that the Einstein upper limit was 1.38 X 10" 13 erg cm" 2 
s" 1 ). Also, fx/hol = 5.35 x 10~ 9 (6 times lower than the Einstein upper limit). Hence a Cas 
is actually a coronal star (as the He I already suggested)! The CDL must be moved to the right 
(slightly) in the bright giant regime. 

3 He I Calculations 

He I A10830 synthetic spectra are generated with the PANDORA stellar atmospheres code both with 
and without the external (i.e., coronal) X-ray field. PANDORA solves the radiative transfer and 
statistical equilibrium equations in a self-consistent, manner with the equivalent 2-level approach 
(see Vernazza, Avrett, & Loeser 1981). Calculations can be made in either plane-parallel or 
splierically-symmetric geometries as either a semi-infinite atmosphere or a finite slab. Macroscopic 
gas velocity fields can be incorporated into the line source function, line profile, and net rate 
equation calculations. External radiation fields incident on the atmosphere can be included in the 
calculations. 

By comparing the various calculations, the influence of a coronal X-ray flux on the He I triplet 
can be explored. We follow the procedure described by Avrett et ai (1976) and Cuntz & Lut- 
termoser (1990) to calculate the excitation and' ionization of helium. The photoionization rate is 
given by 

r o a (u) 

R = 4tt / J u dv , 

Jv i hv 

a(i/) is the photoionization cross section and J u is the mean intensity of both the photosphere/chromo- 
sphere radiation field and the external X-radiation field deduced from the ROSAT observations. 

The following radiative transfer calculations are made in sequence with output from the pre- 
ceeding used as input for the next stage. (1) A 3-level hydrogen model atom is used in a static, 
plane-parallel medium. Ha is treated explicitly and Ly-a and Ly-/3 are assumed to be in detailed 
balance. Balmer and Paschen photo-ionization rates are approximated with an input radiation 
temperature. (2) A 5-level hydrogen atom is used in a static (and where appropriate, dynamic), 
spherically symmetric medium. All transitions are handled in detail. (3) The following species 
will first be calculated assuming a static, plane-parallel medium and second (where appropriate) 
a dynamic, spherically symmetric medium: C I (9-levels), C II (7-levels), O I (13-levels), Mg I 
(7-levels), Mg II (6-levels), Ca I (8-levels), and Ca II (5-levels). (4) A 13-level neutral helium 
atom and a 6-level singly ionized helium atom for both a plane-parallel and spherically symmetric 
medium. All photoionizations and recombinations to the given levels are included. Upon conver- 
gence of the radiative transfer and statistical equilibrium equations with PANDORA, He I A10830 
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synthetic line profiles from semi-empirical models (with and without X-ray fluxes included) are 
compared to the observations of O’Brien and Lambert (1986) and Lambert (1987). These com- 
parisons are a strict test for the X-ray photoionization process in the formation of He I A10830 in 
these cool giants. 

The test of the photoionization method is to produce He I A10830 with W\ % 230 mA for 
a Cas and 20 — 30 mA for a Tau. The incident X-radiation field must not affect the H Balmer 
lines. Must also not affect any other chromospheric emission features (i.e., Mg II h & k and Ca II 
H & K). With the upper limit set by ROSAT for a Tau, assuming an ISM hydrogen column 
density of 8.2 x 10 19 cm 2 and a Raymond— Smith X— ray spectrum, no He I feature formed 
in our representative chromospheric model of a Tau (Kelch et al. 1978). The ROSAT X-ray flux 
measured for a Cas produced a He I A10830 feature at 210 mA(with N H - 2.6 x 10 2 ° cm -2 ), which 
is close to the observed value for this star (see Figure 5)! 

4 Can Stochastic Shocks Produce He I A10830? 

Since the new fx upper limits set by ROSAT for the non-coronal stars are too low to produce 
He I A10830, we need some other process to form this feature. Cuntz & Luttennoser (1990) have 
shown that a 2nd method can give rise to the He I feature. Stochastic shocks can heat the gas 
in the outer regions of these stellar atmospheres to T > 40, 000 K, where sufficient numbers of 
electrons can be pumped to the 2s 3 S metastable state (-20 eV). Cuntz (1987) has shown that the 
dissipation of mechanical wave energy leads to an increase of the mean atmospheric temperatures 
to chromospheric values in the outer layers. Either periodic or non-periodic waves can be modeled, 
where the non-periodic waves change stochastically over a short period. 

A non-periodic wave model of Arcturus (Cuntz 1987) shows a quite interesting behavior. Due 
to the stochastic change of the wave periods, shocks with different strengths are produced. The 
speeds of these shocks are usually different which leads to an overtaking and merging of shocks. 
Since an overtaking shock has a combined strength and therefore an increased speed, it overtakes 
more and more shocks in front of it and attains a very large strength. Such a phenomenon is called 
shock merging. 

The momentum transfer of these strong shocks cause episodic mass loss. Furthermore, the 
energy dissipation of the strong shocks produce high temperatures in the post-shock regions. Oc- 
casionally, T > 40,000 K are reached. Stochastic acoustic wave models for the chromospheres of 
late-type giant stars could lead to a possible explanation for the presence of the He I A 10830 line. 
The highly variable He I profiles in the non— coronal K giants are strong evidence for the presence 
of atmospheric outflows and inflows. The stellar wind flows in the stochastic wave models of Cuntz 
often show the same behavior. 

Cuntz & Luttennoser (1990) have carried out a series of NLTE radiative transfer calculations 
for He I based on hydro models representative of the archetypical, non-coronal star a Boo. We were 
able to produce He I profiles at the proper strength! We had some rather focused requirements: we 
needed shocks at T > 40,000 K that are not too low in the atmosphere (/i % 2 x 10 6 km above 
the visual-continuum formation depths) so that H, Mg II, & Ca II are unaffected. We also needed 
shocks at T > 40,000 K that are not too high (h % 3 X 10 6 km above the visual-continuum 
formation depths) so that densities are high enough to allow sufficient collisional excitations to 
populate the 3 5 state. It was found that radiative transfer effects in a moving medium weaken 
this line with respect to the static calculations. We were unable to produce any He I emission at 
any phase. This lead us to the speculation that long period acoustic waves are also needed in the 
hydrodynamic calculations. 
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for this star. 
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5 Conclusions 


To summarize, this ROSAT program has resulted in the X-ray detection of a previously classified 
non-coronal star ! This detection of a Cas corresponds F x = 0.002F*(O) and is sufficient to form 
the He I A10830 feature in this star. The CDL must be moved to the right on the H-R diagram in 
the bright giant regime! Upper limit to the X-ray flux of a Tau has been reduced by a factor of 70 
over the Einstein value. This flux is insufficient to populate the 3 S state, hence He I feature seen 
in this star cannot arise from the photoionization mechanism. The He I feature in the non-coronal 
stars can be formed by stochastic shocks propagating through their atmospheres. This work has 
been presented in a paper that has just been submitted to the Astrophysical Journal. 
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